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[0001] METHOD FOR ESTIMATING 

SIGNAL MAGNITUDE, NOISE POWER, AND 
SIGNAL-TO-NOISE RATIO OF RECEIVED SIGNALS 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] Tliis application claims priority from U.S. Patent Application Serial No. 

60/501,302, filed September 9, 2003, which is incorporated by reference as if fully set 

forth. 



[0004] FIELD OF INVENTION 

[0005] The present invention relates generally to communication systems and 
techniques which employ quadrature amplitude modulation (QAIM). More specifically, 
the present invention relates to techniques for accurately demodulating an IM-ary QAM 
signal at a receiver. 

[0006] BACKGROUND 

[0007] At present, IVT-ary QAM (M-QAM) is used to provide high-speed data 
transfer in conjunction with a number of mobile data communication applications. 
However, a primary drawback of M-QAM is its computational complexity. One prior 
art approach for addressing this shortcoming involves the use of square M-QAM 
constellations, such as 16QAM, 64QAM, and 256QAM. Square constellations may be 
conceptualized as including two independent q-ary amplitude shift keying ( q -ASK) 
modulation signals as real and imaginary parts (in-phase and quadrature-phase 
components of a modulated carrier), where q = ^f^^ . 

[0008] For example, a 16-QAM signal includes two independent real and 
imaginary 4-ASK signals. A q -ASK signal takes on values selected from the set 
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± — ,±—,...,±- — L Representing a square-constellation M-QAM signal in 9 -ASK 
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format is usually preferred because the complexity of processing a q -ASK signal is 

significantly less than that of processing a general M-QAM signal. 

[0009] A received M-QAM signal (rk) may be represented by the following 

expression: 

r^t = i4*djt +njt ; Equation (1) 

where A is the signal amplitude; rf^ represents a set of data S3niibols encoded onto a 

transmitted M-QAM or q-ASK signal; and represents noise. While for M-QAM 

signals, and are complex quanities, forg-ASKsignals they are both purely real. 

[0010] In order to demodulate M-QAM and q-ASK signals, it is necessary to 

determine the values of one or more signal parameters at the receiver, such as 
amplitude, noise power, and signal-to-noise ratio (SNR). Even in the case of a (y-ASK 
signal representation, which provides enhanced computational simplicity relative to M- 
QAM representation, existing techniques for estimating these signal parameters are 
inefficient and computationally complex. 

[00 1 1] An improved technique for estimating signal parameters at the receiver is 

needed to provide efficient demodulation of M-QAM and g-ASK signals. 

[0012] SUMIMARY 

[0013] The present invention provides an improved system and method for 
estimating one or more parameters of a received M-QAM or q-ASK signal, such as 
amplitude and SNR. A first embodiment of the invention estimates the amplitude of an 
M-QAM signal based upon known or ascertainable phase information regarding a 
plurality of transmitted symbols. A respective set of received sjnnbols, corresponding 
to the plurality of transmitted symbols, is recovered. Each of the plurality of received 
symbols is multiplied by a complex unit vector with a phase that is opposite in sign to 
the complex transmitted data symbol to generate a set of products. The set of products 
is summed, and the real part of the sum of products is then determined. The absolute 
values of the Imown transmitted sjnnbols are summed to generate a total magnitude 
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value. The real part of the sum of products is divided by the sum of transmitted 
magnitude values to generate an estimate of the amplitude of the M-QAM signal. 
[0014] A second embodiment of the invention estimates the amplitude of a q-ASK 
signal at a receiver based upon known or ascertainable magnitude information 
regarding a plurality of transmitted symbols. A respective set of N received samples 
corresponding to the transmitted symbols is recovered. For each of the N samples, the 
sample is multiplied by the corresponding sign of the transmitted sjnnbol to generate a 
set of products. This set of products is then summed to generate a first sum. The 
absolute values of the transmitted symbols are then summed to generate a second sum. 
The first sum is divided by the second sum to generate an estimated amplitude for the 
q-ASK signal. 

[0015] A third embodiment of the invention estimates the amplitude of an M- 
QAM or q-ASK signal with no assumed knowledge of the specific transmitted data 
symbols. Only statistical Imowledge of the amplitude of the transmitted data symbols 
is assumed. A respective set of received samples corresponding to the transmitted 
symbols is recovered. The absolute values of the received samples are determined and 
summed to generate a first sum. The mean of the absolute values of the amplitudes of 
transmitted S3nnbols is determined and then multiplied by the total number of received 
symbols to generate a product. The first sum is divided by the product to generate an 
estimated amplitude for the M-QAM or q-ASK signal. 

[0016] Other embodiments of the present invention utilize second-order and 

fourth-order moments of received samples, a maximum likelihood searching process, or 
a Kurtosis estimation process to estimate amplitude, noise power, and signal-to-noise 
ratio of a received signal. 

[0017] BRIEF DESCRIPTION OF THE DRAWING(S) 

[0018] FIG. 1 is a flow diagram of a procedure for estimating the amplitude of an 

M-QAM signal based upon Imown or ascertainable phase information concerning a 
plurality of transmitted sjonbols, according to a first embodiment of the invention. 
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[0019] FIG. 2 is a flow diagram of a process for estimating the amplitude of a q- 

ASK signal at a receiver based upon known or ascertainable magnitude information 
regarding a plurality of transmitted S3niibols, pursuant to a second embodiment of the 
invention. 

[0020] FIG. 3 is a flow diagram of a procedure for using absolute values of 
received samples to estimate the amplitude of an M-QAM or q-ASK signal, pursuant to 
a third embodiment of the invention. 

[0021] FIG. 4 is a flow diagram of a procedure for using second-order and fourth- 
order moments of received samples to estimate amplitude, noise power, and signal-to- 
noise ratio of an M-QAM signal that includes a set of transmitted symbols, pursuant to 
a fourth embodiment of the invention. 

[0022] FIG. 5 is a flow diagram of a procedure for using second-order and fourth- 
order moments of received samples to estimate amplitude, noise power, and signal-to- 
noise ratio of a q-ASK signal that includes a set of transmitted symbols, pursuant to a 
fifth embodiment of the invention. 

[0023] FIGS. 6A, 6B and 6C are a flow diagram of a procedure for using a 
maximum likelihood searching process to estimate the amplitude, noise power, and 
signal-to-noise ratio of a q-ASK or IM-QAM signal that includes a set of transmitted 
s3nnbols, pursuant to a sixth embodiment of the invention. 

[0024] FIGS. 7A and 7B are a flow diagram of a Kurtosis estimation process for 

estimating the signal-to-interference ratio of an M-QAM or q-ASK signal using second- 
order and fourth-order moments of received samples, pursuant to a seventh 
embodiment of the invention. 

[0025] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 

[0026] Improved systems and methods are disclosed for estimating one or more 

parameters of a received M-QAM or q-ASK signal, such as amplitude and signal-to- 
noise ratio. A first embodiment of the invention is a method 100 which estimates the 
amplitude of an M-QAM signal based upon Imown or ascertainable phase information 

concerning a plurality of transmitted s3nnbols (dk). If the receiver has knowledge of 
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the transmitted M-QAM S3nnbols , for example in systems employing pilot signals, 
this information may be used to estimate the amplitude, A, as: 

Reff;r,exp[-yW*)]] 
A = , ^ ; Equation (2) 

where e{df^ ) is the phase of complex symbol and N is the number of received symbols 

being processed. Re( ) represents the Real operator, noting that one may substitute 
many other operators without substantively changing the invention. 
[0027] More specifically, referring to FIG. 1, a respective set of received symbols 

(rk), corresponding to the plurality of transmitted symbols, is recovered (step 101). 
Each of the plurality of received sjrmbols (rk) is multiplied by exp [-je(dk)] to generate a 
set of products (step 103), where 0(dk) represents the phase of a corresponding 
transmitted symbol (dk). The set of products is summed (step 105), and the real part of 
the sum of products is then determined (step 107). The absolute values of the 
transmitted sjnnbols | (dk) | are summed to generate a magnitude value (step 109). The 
real part of the sum of products is divided by the magnitude value to generate an 
estimated amplitude for the M-QAM signal (step 111). 

[0028] A second embodiment of the invention is a method 200 which estimates 
the amplitude of a q-ASK signal at a receiver based upon known or ascertainable 
magnitude information regarding a plurality of transmitted symbols (dk). If d^^ is q- 

ASK signal, then amplitude estimation may proceed in accordance with the following 
equation: 

N 

A = . Equation (3) 

With reference to FIG. 2, a respective set of N received samples (yk) corresponding to 
the transmitted s3rmbols (dk) is recovered (step 201). For each of the N samples, the 
sample (yk) is multiplied by a corresponding sign (dk) to generate a set of products 
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(yk)*sign(dk) (step 203). This set of products is then summed to generate a first sum 
(step 205). The absolute values of the transmitted sjnnbols |(dk)| are then summed to 
generate a second sum (step 207). The first sum is divided by the second sum to 
generate an estimated amplitude for the q-ASK signal (step 209). 
[0029] A third embodiment of the invention is a method 300 which uses absolute 

values to estimate the amplitude of an M-QAM or q-ASK signal which includes a set of 
N transmitted sjmabols (dk), where N is a positive integer greater than one. In many 
cases, the receiver lacks information about the transmitted S3nnbols d^^ , such as the 

phases and magnitudes of the symbols. Accordingly, the amplitudes of the received 
samples may be utilized to generate an amplitude estimate for both M-QAM and q- 
ASK signals from the relationship: 

A-^f^y Equation (4) 

where E\d^\ is the mean of the amplitude of signal dj^ and may be computed 
theoretically or empirically. For example, for 4-ASK modulation with e |±^,±^| 
with equal probability, £1^^1 = 1, the amplitude estimate is: 

^^^Zhl- Equation (5) 

[0030] With reference to FIG. 3, a respective set of received samples (yk) 
corresponding to the transmitted symbols (dk) is recovered (step 301). The absolute 
values of the received samples | (yic) | are determined (step 303) and sunmied (step 305) 
to generate a first sum. The mean of the absolute values of the amplitudes of 
transmitted symbols, E | (dk) | , is determined (step 307) and then multiplied by N to 
generate a product, N* E | (dk) | (step 309). The first sum is divided by the product to 
generate an estimated amplitude for the IM-QAIM or q-ASK signal (step 311). 
[0031] Through the use of second-order and fourth-order moments of received 

samples, a method 400 of a fourth embodiment of the invention estimates the 
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amplitude, noise power, and signal-to-noise ratio of an M-QAM signal that includes a 
set of transmitted sjnnbols (dk). The noise at the receiver takes the form, = + jnl , 
where and can be approximated as two independent Gaussian variables with 
zero means and same variance cr,^ E\ni^f =2<7^ which is often exact, or at least a good 
approximation. The amplitude is estimated as: 



i 



2E 



The noise power is estimated as: 



Equation (6) 



&'„=E\rf^A'E\dA\ 



Equation (7) 

Using the above two Equations (6) and (7), the SNR is estimated as: 

SNR = A^E\df^\^ la] . Equation (8) 



[0032] Referring now to FIG. 4, the method 400 begins when a respective set of 
received samples (rk) corresponding to the transmitted S3rmbols (dk) is recovered (step 
401). The mean of the absolute values of the amplitudes of the transmitted s3mibols, 
E I (dk) I , as well as the mean of the absolute values of the amplitudes of the received 
samples, E|(rk)|, are determined (step 403). The estimated amplitude A is 
determined using Equation 6 (step 405). Estimated noise power o\ is determined from 
estimated amplitude A using equation 7 (step 407). Optionally, estimated SNR is then 
determined using Equation 8 (step 409)). 

[0033] Through the use of second-order and fourth-order moments of received 

samples, a method 500 of a fifth embodiment of the invention estimates amplitude, 
noise power, and SNR of a q-ASK signal that includes a set of transmitted symbols (dk). 

When a real q-ASK signal is transmitted, the second and fourth order moments may 
be used to estimate the signal amplitude and power as: 



\3[E{d,'))-E{d,') ' 



Equation (9) 
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^ 3[E(r, )f-£(r, ) ^ Equation (10) 

p[E{d,')[-Eid,') 

The noise power is estimated as: 

= E{r^^)-A^E(di,^) . Equation (11) 

Using Equations (10) and (11), the SNR is estimated as: 

SNR = A^E{d^'')lal . Equation (12) 

[0034] Expectations involving dk may be computed empirically or theoretically. 
For example, for 4-ASK modulation E(d^^) = 1.25 and E(flf/) = 2.5625 . The signal 
amplitude and power may be estimated by substitution with the determined values 



e.g 



^^4 p(g(^.'))'-^(r/) andA^ ^ h(Eir,')f-Eir,') actively. The SNR is estimated 
¥ 2.125 V 2.125 ' ^ 



by using Equation (12). Various methods for generating estimates of the expectations 
involving rk may be used. 

[0035] With reference to FIG. 5, the method 500 begins when a respective set of 

received samples (rk) corresponding to the transmitted symbols (dk) is recovered (step 
501). The mean of the amplitudes of the transmitted symbols, E(dk), as well as the 
mean of the amplitudes of the received samples, E(rk), are determined (step 503). 
Estimated amplitude A is determined using Equation 9 (step 505). Optionally, 
estimated power is determined using Equation 10 (step 507). Estimated noise power 
o^n is determined from estimated amplitude A using Equation 11 (step 509). 
Optionally, estimated SNR may then determined using Equation 12 (step 511). 
[0036] Using a maximum likelihood searching process, a method 600 of a sixth 
embodiment of the invention estimates at least one of the amplitude, noise power, and 
signal-to-noise ratio of a q-ASK or M-QAM signal that includes a set of transmitted 
s3rmbols (dk). Assuming the received samples are sufficiently independent, the joint 
probability density function of ri,r2,.,.,r^ may be assumed to be the product of the 
marginal distributions of rk, that is: 
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N 

f{A, a^\r^, ) = fl /(^' i 'i ) 5 Equation (13) 

where fiA,(T„;r^) is an estimate of the density of rk. The goal of the maximum- 
likelihood amplitude estimator is to find the amplitude, A, that maximizes the joint 
probability density function as shown by Equation (13), i.e.: 

i4 = argmax /(i4,cr„;rpr2,...,r^) . Equation (14) 

A 

[0037] The following search algorithm finds the maximum likelihood estimate for 

amplitude, and optionally for noise power and SNR: First, the received vector is 

1 ^ 

normalized such that its average power is unity, (i.e. — ZlrJ^ = l .) A^^in and A^^^ are 

N k=i 

selected according to the minimum and maximum SNR over which to search as follows: 



SNR^,„ ; and Equation (15) 



^max = j . Equation (16) 



max 



[0038] The amplitude search step size A is selected such that it achieves the 

desired amplitude resolution. Next, A^=A^^-^kA is determined, where 0<k<K 

and K = — . For each amplitude Aj^ , a corresponding noise power is determined 

A 

as: 

cr„^ =1-£|^^|'a/ . Equation (17) 

[0039] The joint probability density function /(i4^,cr„;rpr2....,r^) is calculated for 
each Af^. The value /i^ that corresponds to the maximum joint PDF, (i.e., 
A= argmax /(i4t,(7„;ri,r2,...,r^) ), is searched for, and then A is the estimated amplitude. 

[0040] Additionally, (and optionally,) a value for the estimated noise power can 
be determined from the estimated amplitude A as the total power less the product of 
the amplitude estimate squared and the expected data sjnnbol power, (i.e., o\ = E(rk)2 
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- A2 E(dk)2). Optionally, a value for the estimated SNR can be determined by taking 
the ratio of the estimated signal power and the estimated noise power, (i.e., from the 
relationship SNR = [A2 * E(dk)2] / aV) 

[0041] With reference to FIGS. 6A, 6B and 6C, the method 600 briefly described 
in the foregoing paragraphs is described in greater detail. A respective set of received 
samples (rk) corresponding to the transmitted symbols (dk) is recovered (step 601). N 
received samples (rk) are normalized, (N being a positive integer greater than one), 

such that — Z|rJ^ =1 (step 603). Trial values for A^^ and A^^^ are selected according to 

N k^i 

the minimum and maximum SNRs over which to search using Equation 15 and 16 
(step 605). 

[0042] An amplitude search step size A is selected to achieve a desired or 

predetermined amplitude resolution (step 607). Next, a value for Aj^ = +/:A is 

determined, where 0<k<K and K = --^ ^ (step 609). For each amplitude A, , a 

A 

corresponding noise power is then estimated using Equation 17 (step 611). 

[0043] A joint probability density function /(A^ , a„ ; ) is then determined 

for each (step 613). 

[0044] For both M-QAIVI and q-ASK signals, a main objective of the maximum- 
likelihood searching process is to find a value for the estimated amplitude, A , that 
maximizes the joint probability density function using Equation 14 (step 615). At step 
617, a search is conducted to find a value for A^ that corresponds to the maximum joint 
probability density function. This value for A^ is substituted into the Equation 
A= argmax /(A4,cr„;rpr2» to obtain a value for A, representing an estimated 

amplitude value (step 619). 

[0045] Optionally, a value for estimated noise power can be determined from the 

estimated amplitude A as the estimated total power less the product of the amplitude 
estimate squared and the expected data symbol power, (i.e., o^n = E(rk)2 - E{dk)^ 
(step 621)). Additionally, (and optionally), a value for estimated SNR can be 
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determined as a ratio of the estimated signal power and the estimated noise power, 
(i.e., SNR = [A2 * E(dk)2] / (step 623). 

[0046] Pursuant to a method 700 of a seventh embodiment of the invention, a 
Kurtosis estimation process estimates the signal-to-noise ratio of an M-QAM or q-ASK 
signal using second-order and fourth-order moments of received samples (rk). As a first 
step, the Kurtosis of the noise only (Kcg) and signal only (Ksig) present is computed. In 
this example, complex Gaussian noise is used for illustrative purposes. However, other 
noise distributions may be assumed and do not deviate substantively from this 
invention. Expectations involving dk and Uk may be computed empirically or 
theoretically, e.g., the Kurtosis values are determined as follows: 



^ ^IKI I = 2 ; and Equation (18) 

K,^-E^^; Equation (19) 

where complex, circularly symmetric Gaussian noise is assumed. For example, with 
QPSK data, Ksig = Kqpsk = 1, for square 16QAM data, Ksig = Kigqam = 1.32. The Ksig 
parameter for any signal constellation can be readily obtained. 
[0047] The Kurtosis for rectangular and 'cross' QAM constellations can also be 
used with this method and determined as outlined above. The Kurtosis may be written 
in terms of, and solved for the SNR as follows: 

Any of a variety of Kurtosis estimates or approximations may be used. 
[0048] Referring to FIGS. 7A and 7B, a flow diagram of the details of the 
foregoing Kurtosis process is shown. As a preUminary matter, the second-order 
moment is defined as£{|/ip}= £{|ni|^}+£{|flft|'}, and the fourth-order moment is defined 
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as E' 



Ik* r l'^ ^(K r K ^(H* r K ^^Ih* r I^IM* r where dk denotes the transmitted symbols and 



nk denotes a noise component that is recovered with the received samples rk (step 701). 
At step 703, the fourth-order moment is divided by the second-order moment so as to 
implement a Kurtosis operation as follows: 



attributable to noise alone (step 705). For example, assuming the complex circularly 
symmetric Gaussian noise, the Kurtosis of the noise alone is determined using 
Equation 18. The Kurtosis attributable to the signal alone, (Ksig)^ is determined using 
Equation 19 (step 707). For example, with QPSK data, Ksig = Kqpsk = 1 may be used. 
For square 16QAM data, Ksig = Kigqam = 1.32 may be used. For high order modulation, 
square M-QAM, the signal Kurtosis Ksig approaches that of the square complex uniform 
distribution: Hm {Kj^_qj^ ) = 1.4 . However, use of a square constellation is not required, 

since the Kurtosis for rectangular and 'cross' QAM constellations can also be utilized in 
conjunction with this method and pre-computed in the same way. 
[0050] Once an appropriate value for Ksig (e.g., Km-qam in the case of a M-QAM 

signal) is determined, the method 700 of FIG. 7 advances to step 709 where the signal- 
to-noise ratio (SNR) is estimated using Equation 20. 




Equation (21) 



[0049] 



In order to estimate the SNR, it is necessary to determine the Kurtosis 



* 
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